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Abstract 

We calculate several correlation functions and distribution functions of dynamical 
variables for the gauge glass and the Villain model using the spin wave approxima- 
tion and the gauge transformation. The results show that the spin wave approxi- 
mation gives the exact solutions on the Nishimori line in the phase diagram. This 
implies that vortices play no role in the thermodynamic behavior of the system as 
long as some correlation and distribution functions are concerned. These results 
apply to any dimensions including the two-dimensional case. 
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1 Introduction 

The gauge glass, which is an extension of the XY model to the random spin 
system, attracts both theoretical and practical interest. This model provides 
good descriptions of some physical situations, for example, the XY magnet 
with random Dzyaloshinskii-Moriya interactions [1] and a granular supercon- 
ductivity and Josephson-j unction arrays with positional disorder [2] [3]. Espe- 
cially, the two-dimensional gauge glass has been studied actively because of the 
interesting relation between the Kosterlitz-Thouless (KT) transition and the 
effect of disorder. This relation was first discussed by Rubinstein, Shraiman 
and Nelson [1] . They have shown that for small amount of randomness, as the 
temperature is decreased, there appears first the paramagnetic phase, then 
a KT-like phase, and finally again a paramagnetic phase. Thus a reentrant 
transition appears in the gauge glass. After their work, Natterman, Scheidl, 
Korshunov and Li [4] corrected their description of the reentrant transition 
to the one without reentrance. These latter authors pointed out an overesti- 
mation of vortex pair density in the previous work. Furthermore, Sheidl [5] 
found a new ordered phase in which single vortex excitations occur in the 
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reentrant phase, and Maucourt and Grempel [6] suggested from Monte Carlo 
simulations that there is no indication of a low-temperature reentrant phase. 
Today, the two-dimensional gauge glass at low temperatures is generally be- 
lieved to have no reentrance [7], except for some researchers [8]. Hence, the 
present major interest in the two-dimensional gauge glass is the structure of 
the KT phase, namely, the boundary between the freezing phase of the vortex- 
pair excitations and the non-freezing phase. 

The gauge glass has also been investigated in the field of the spin glass theory 
because the gauge glass is an extension of the Ising spin glass with continuous 
spin variables [9] . Particularly, the method of gauge transformation, which was 
developed in the study of the spin glass, is a powerful technique for deriving 
analytical results for gauge glass. Ozeki and Nishimori [10] found the exact 
solution of the internal energy of the gauge glass under a special condition 
using the method of gauge transformation. This special condition relevant to 
the exact solution corresponds to a line in the phase diagram called the Nishi- 
mori line. Ozeki and Nishimori also showed that the phase boundary between 
the ordered and the disordered phases runs parallel to the temperature axis. 
Although they could obtain the exact solution on the Nishimori line, there 
exists a mysterious property associated with the line: For example, the exact 
solution of the internal energy on the Nishimori line has no singularity as a 
function of temperature although the the line runs across the phase boundary. 
From these facts in mind, we aim to clarify what occurs under the special 
condition of this line. 

In this paper, we calculate some gauge invariant quantities on the Nishimori 
line. First, we define the gauge glass and show that the exact solution of 
the internal energy can be calculated under the special condition in Section 
2. Next, we introduce the gauge invariant correlation functions and calculate 
them using the spin wave approximation, gauge transformation and the Villain 
model in Section 3. The gauge invariant distribution functions are introduced 
and calculated in Section 4. From the results of Sections 3 and 4, we discuss 
the physical properties of the Nishimori line in Section 5. 



2 General properties 

2.1 Gauge glass 

The gauge glass is defined by the Hamiltonian 
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where Aij is the quenched random phase shift and the couphng constant J is 
positive. The sum runs over all nearest-neighbor pairs on a lattice. Note that 
we do not specify the spatial dimensionality explicitly here. We assume that 
Aij is independently distributed at each bond as follows, 



where the magnitude of Kp controls the tendency toward ferromagnetism; the 
system becomes the pure ferromagnetic XY model in the limit of large Kp. 
The function Io{x) appearing in the denominator of eqn (2) is the modified 
Bessel function /o(a;) = /g^'^d^ ^xcose^ Note that we express the thermal aver- 
age by angular brackets (■ ■ ■) and the average over the distribution of quenched 
randomness by square brackets [■■■]. 



2.2 Exact solution of the internal energy 



The gauge transformation is a powerful method to calculate the exact solution 
of the internal energy for the gauge glass [11]. We introduce the gauge trans- 
formation to our model and derive the exact internal energy under a special 
condition. 

The gauge transformation is defined as 



9i^6i-(f)i, Aij ^ Aij - 4)i + (j)j, (3) 

where 0j is a gauge variable fixed arbitrarily at each site. The Hamiltonian (1) 
is invariant under this transformation. The internal energy is written explicitly 

as 



I- A A ■ P^V ^1-] I 'ZTT 



cos(6»i-6»j-Aij) 



cos(ei-ej-Aij) 



where the ranges of integration over Aij and 9i are from to 27r. Here we 
apply the gauge transformation to eqn (4). Since the gauge transformation is 
just a change of running variables, the value of eqn (4) is independent of the 
choice of {(^j} and therefore we may integrate eqn (4) over all gauge variables 
each from to 27r, and divide each integration of (^j by 27r. Using the 
special condition K = Kp, we obtain the exact internal energy as 
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Fig. 1. The schematic phase diagram and the Nishimori Une (dashed). The spin 
glass phase exists in three dimensions or over. 



Here the total number of interacting pairs is denoted as Nb- The function 
Ii{x) is the modified Bessel function = /Q^'^d6'cos^ gxcose _ 
We could carry out the calculation of the internal energy exactly owing to 
the gauge invariance of the Hamiltonian and the special condition K = Kp. 
The gauge invariance plays an important role in the gauge glass or random 
spin system because, for the gauge invariant quantities, the exact solution 
similar to the case of the internal energy can be calculated by using gauge 
transformation. The condition K = Kp defines a line (the Nishimori line) in 
the phase diagram as in Figure 1. Although the line K = Kp clearly runs 
across the phase boundary, the exact solution of the internal energy (5) is a 
simple function because the modified Bessel functions has no singularity. Thus, 
we may consider that the singular part of the internal energy vanishes on the 
Nishimori line. The same behavior of the internal energy as in the gauge glass 
appears in the Gaussian spin glass model and ±J spin glass model [11]. 
We can further evaluate some quantities using the spin wave approximation, 
which is a characteristic method in continuous spin systems, and obtain some 
information not available in Ising spin systems, which helps us to understand 
the significance of the Nishimori line. In particular, it is believed that the 
vortices created by the periodicity of eqn (1) cause the singularity relevant 
to the phase transition in the two-dimensional gauge glass. Accordingly, we 
expect that a new physical picture relevant to the Nishimori line may be 
found by using the characteristic method in the continuous spin system. This 
program is carried out in the next section. 




(5) 
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Fig. 2. Schematic diagram of the correlation function Cp{r) and chirahty. Left figure 
shows a path from site to site r. The right figure shows two states with different 
chiralities. 

3 Gauge invariant correlation functions 



3. 1 Definitions 



In order to investigate the physical significance of the result in the previous 
section, we introduce the gauge invariant correlation functions as 



Cle - 
Cch ■ 



exp i{0Q -Or-Y^ Ak-i,k) 

\ k=l 
C0s(6'o — 9l — Aoi) COs{9r — Or+l — 

Sm{6o - 6*1 - Aoi) Sm{9r - 6r+l - ^r,r+l) 



(6) 
(7) 



Equation (6) is the gauge-invariant correlation function which measures the 
correlation along a path from site to site r taking the phase twist of Aij into 
consideration. In Figure 2, we sketch a schematic picture of a path of corre- 
lation function Cp. Equation (7) means the local energy correlation function. 
The local energy correlation function is also gauge invariant. Equation (8) rep- 
resents the chirality correlation function. Chirality introduced by Villain [12] 
is an Ising-like degree of freedom. The chirality distinguishes the two states 
which are transformed to each other using the mirror transformation (Figure 
2). The chirality correlation function is also gauge invariant. 
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3.2 Spin wave approximation 



Let us now consider the Hamiltonian introduced in eqn (1) in two dimensions. 
At sufficiently low temperatures, the effect of thermal fluctuations is small, so 
that the spins tend to align to each other with the twist Aij taken into account. 
Hence, we are able to expand the cosine term near {9i — 9j — A^j) = 0, 



Hc^-Y.iO^-0,-A,r (9) 

We note that this approximation neglects periodicity of cosine in eqn (1). It 
is generally believed that, if we want to describe the phase transition in two- 
dimensional gauge glass, we must consider the contribution of periodicity or 
vortices [7]. To ensure the validity of the expansion of eqn (9), we must assume 
that the phase twist A^j is sufficiently small, so that eqn (2) is expressed as 



P{Aij) oc e-^-^^^/^ (10) 

The constant is equal to the variance of the Gaussian distribution [A'^j] = 
a. In the same way as in the spin wave approximation, the distribution of eqn 
(10) neglects periodicity of Aij. This approximation of the distribution of Aij 
is expected to be sufficient to describe correlations at low temperature. 

Using the spin wave approximation of eqns (9) and (10), we can calculate the 
gauge invariant correlation functions as 



where Cor and Gq^^ are the lattice Green functions. 



_ly^ l-cosfc-r 

~ iV V 4 - 2 cos - 2 cos K 



1 


1 


K'p 


~ K 


1 


1 




~ H 



q{±) _ 1 (1 — COS A;i)(l ± cosA; ■ r) 



(15) 
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Now, applying the special condition K = Kp, we can obtain these correlation 
functions as follows, 

Cp = e-^"'-/2, Cie = e^", Cch = 0. (16) 

Surprisingly, the chirality correlation function vanishes under the condition 
K = Kp. The chirality degrees of freedom behave perfectly independently 
from place to place. 



3.3 Gauge transformation 

From comparison of the results of the spin wave approximation with the exact 
results, we can investigate the difference between the model with periodic- 
ity and the model without periodicity. In this section, we calculate the exact 
solutions of Cp, Cie and Cch for models with periodicity using gauge trans- 
formation under the special condition K = Kp. After the derivations of the 
exact solutions of correlation functions, we study the asymptotic forms of the 
exact solutions at low temperatures and compare them with the solutions in 
the previous section using the spin wave approximation. 
The correlation functions which we introduced in Section 3.1 are invariant un- 
der the gauge transformation, and therefore the same calculation method as 
for the internal energy (5) in Section 2 can be applied to the calculation of the 
correlation functions. As a consequence, we obtain the following expression 
under the special condition K = Kp, 

When K = Kp, the chirality correlation vanishes exactly and this solution 
is consistent with the result of spin wave approximation (16). At low tem- 
perature, the asymptotic behaviors of Cp and Cie in eqn (17) are written as 
follows, 

Cp~e-^"^/^ Cie-e-^'^'-". (18) 

Equation (18) is also consistent with the solution of spin wave approxima- 
tion expressed in eqn (16). This comparison suggests that under the special 
condition K = Kp, the system forms a characteristic structure which is not 
influenced by the vortices. The analysis of the Villain model in the following 
reinforces this picture. 



C 



ch 



0. 



(17) 
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3-4 Villain model 

The Hamiltonian of the spin wave approximation (9) neglects periodicity in- 
cluded in the cosine term (1). If we need to describe the system at low tem- 
peratures more accurately, it should be necessary to consider periodicity in 
the Hamiltonian. Accordingly, we consider the Villain model which adds the 
periodicity to the spin wave Hamiltonian 



where rriij is an integer between — oo and oo. 

The random phase shift (10) also neglects the periodicity of modulo 27r. There- 
fore we also add 2mij7i into the probability distribution of Aj (10). The con- 
figurational average is described as 



From eqns (19) and (20) and the gauge transformation, we can calculate the 
correlation functions discussed in the previous section. We first consider the 
gauge-invariant correlation function Cp. Applying the gauge transformation 
and the condition K = Kp, we obtain the following form. 




(19) 




(20) 




(21) 







To integrate over A^j, we use the following relation. 




(22) 
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Substituting eqn (22) into eqn (21), we obtain the gauge-invariant correlation 
function as 



C7p = e-^^"'^/^ (23) 

where we have used the relation K = Kp. The result of this calculation agrees 
with the one of the spin wave approximation (16) and the asymptotic form 
of the exact solution (18) under the condition K = Kp. The local energy 
correlation function and the chirality correlation function are derived by the 
same calculation as in the gauge invariant correlation. These results are written 
by simple formulas as 



ae = e-^", Cch = 0. (24) 

Those two correlation functions also agree with the results of the spin wave 
approximation (16) and asymptotic exact solution (18) under the condition 
K = Kp. We may conclude from these results as follows. The spin wave 
approximation neglects periodicity, and the results of the spin wave approxi- 
mation do not include the effects of vortices. The results for the Villain model 
in contrast include the effects of vortices. Accordingly, the agreement of the 
results of spin wave approximation with those of the Villain model suggests 
the absence of the effects of vortices on the correlation functions we calculated 
under the condition K = Kp. Vortices are perfectly irrelevant to the gauge 
invariant correlation functions. 



4 Distribution functions 



Irrelevance of vortices in the gauge invariant correlation functions on the Nishi- 
mori line can be verified by using a distribution function with gauge invari- 
ance. This gauge invariant distribution function measures the distribution of 
the random phase twist between site i and j. The definition is 



Fix) 



(25) 



Using the spin wave approximation and the Villain model, we obtain the 
comparable results as in the gauge invariant correlation functions when K = 
Kp. The results using the spin wave approximation and the Villain model are, 
respectively. 
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Psw(a:) = ^§exp(-^x2) (26) 

^vmain(x) = ^exp(-^x2), (27) 

the latter being the exact solution. These two expressions are completely iden- 
tical. 

The distribution function of two bonds can also be calculated by the spin 
wave approximation and the Villain model. We first show the definition of the 
two-bond distribution function 



P{x,y) 



{S (X - (Ok -01- A,i)) .6iy-iem-0n~ Amn)) ) 



(2J 



Equation (28) is calculated by using the spin wave approximation and the 
Villain model for K = Kp as, again the latter being the exact solution, 



A/iiiain(a;, y) = 




Similarly to eqns (26) and (27), these results agree with each other. We there- 
fore see immediately that vortices have no effects on the distribution function 
(28). It also follows that two-bond variables behave independently because 
eqns (29) /(30) are just products of eqns (26)/ (27). 
It is instructive to consider another two-bond extension. 



P'{x,y) 



{6 {x - {9k -01- Am)) )-{S{y- {d^ - - A^„,)) ) 



(31) 



Although eqn (31) is gauge invariant, we can not calculate the exact solution 
even if we apply gauge transformation to eqn (31). Nevertheless, we can obtain 
the result by the spin wave approximation as 



where A = Gkn — Gim — Gkm + Gin- In contrast to eqn (28), the distribution 
function (31) can not be written as the products of eqn (26) under the con- 
dition K = Kp. Bond variables are correlated in this sense. Inequivalence of 
eqns (28) and (31) is a highly non-trivial fact. If we measure the two-bond 
correlation within a given system as in eqn (28), the result becomes uncorre- 
lated: P{x,y) = P{x)P{y). If we, on the other hand, look at the correlation 
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of the thermal averages of single-bond distribution function at two different 
locations as in eqn (31), these are correlated: P'{x,y) ^ P{x)P{y). Physical 
significance of this counter-intuitive results needs further clarification. 



5 Summary and discussions 

In this paper, we have investigated the properties of the gauge glass on the 
Nishimori line through the gauge invariant quantities. Firstly we calculated 
the gauge invariant correlation functions using the spin wave approximation, 
gauge transformation and the Villain model. From these calculations, we found 
that the results of the spin wave approximation are identical to the results by 
the method taking periodicity into account. Thus, we expect that the peri- 
odicity or vortices has no influence on the system (gauge invariant quantities 
to be precise) under the condition K = Kp. The effect of vortices can be ne- 
glected on the Nishimori line. 

Next, we calculated the gauge invariant distribution functions using the spin 
wave approximation and the Villain model. We found the same results as for 
the gauge invariant correlation functions. Moreover, we showed that the dis- 
tribution function (28) is independently distributed at each bond. An variant 
of eqn (28) was also calculated by using the spin wave approximation. 
We may conclude from these results as follows. In the two-dimensional gauge 
glass, the effect of vortices plays an important role because the phase transi- 
tion of the model is described by vortex-pair unbinding. Accordingly, the spin 
wave approximation which neglects the effect of vortices is considered in gen- 
eral not to be able to describe the singularity related to the phase transition. 
However, our analysis based on the gauge invariant quantities predicts that 
under the condition K = Kp, the singularity related to the phase transition 
vanishes for gauge invariant correlation functions even if we consider the effect 
of vortices in the system. Thus, the renormalization group arguments [7] [13], 
which analyze the effects of vortices asymptotically, should be reconsidered 
seriously under the present perspective. Further investigation of the property 
of the Nishimori line is necessary. 
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